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I.  INTRODUCTION 


A  major  area  of  concern  in  shell  design  is  the  total  aerodynamic  drag.. 
The  designer,  ever  desirous  of  increasing  the  range  and/or  terminal  velocity 
of  projectiles,  is  eager  to  decrease  the  aerodynamic  drag. 

The  total  drag  of  projectiles  can  be  divided  into  three  components:  (1) 
pressure  drag  {excluding  the  base  region),  (2)  viscous  (skin  friction)  drag, 
and  (3)  base  drag.  For  a  typical  shell  at  M  =  .90  the  relative  magnitudes  of 
the  aerodynamic  drag  components  are:  (1)  pressure  drag,  20%,  (2)  viscous 
drag,  30%,  and  (3)  base  drag,  50%.  The  pressure  and  viscous  components  gener¬ 
ally  cannot  be  reduced  significantly  without  adversely  affecting  the  stability 
of  shell.  Recent  attempts  to  reduce  the  total  drag  are  therefore  directed  at 
reducing  the  base  drag. 

A  number  of  studies  have  been  made  to  examine  the  total  drag  reduction 
due  to  the  addition  of  a  boattail.  Although  this  is  very  effective  in  reduc¬ 
ing  the  total  drag,  it  has  a  negative  impact  on  the  aerodynamic  stability  of 
shell  especially  at  transonic  velocities.  An  excellent  review  of  the  effect 
of  boattailing  on  total  drag  and  base  pressure  is  presented  in  Reference  l. 

Another  effective  means  of  reducing  the  base  drag  is  that  of  'base  bleed' 
or  'base  injection.'  In  this  method,  a  small  amount  of  mass  is  injected  into 
the  base  region  which  increases  the  base  pressure  and  thus  reduces  the  base 
drag.  Recent  range  and  precision  tests2  of  a  155mm  projectile  with  and  with¬ 
out  base  bleed  have  been  conducted  and  an  85%  reduction  in  base  drag  was 
obtained.  Presently  the  XM864  is  an  active  projectile  design  which  is  attempt¬ 
ing  to  use  the  base  bleed  concept  for  increased  range.  This  concept  of  mass 
injection  at  the  projectile  base  has  been  widely  studied  for  supersonic  flows 
and  much  of  the  work  has  been  reported  in  Reference  3.  One  limited  study  at 
supersonic  speeds  was  made  at  BRL  and  the  results  were  reported  by  Dickinson.4 
The  supersonic  regime  has  typically  been  the  area  where  increased  range  due  to 


1.  Sedney,  R.,  "Review  of  Base  Drag,"  Report  No.  1337,  U.S.  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland ,  21005,  October 
1966  (AD  808767). 

2.  "155mm  ERFB  Base  Bleed  Range  and  Precision  Tests,"  Conducted  at  Proof  and 
Experimental  Test  Establishment ,  Nicolet,  Quebec, for  Space  Research 
Corporation,  January  11,  1978 . 

3.  Murthy,  S.N.B.  (Ed.),  "Progress  in  Astronautics  and  Aeronautics:  Aerody¬ 
namics  of  Base  Combustion,"  Vol .  40,  AIAA,  New  York,  1976. 

4.  Dickinson ,  E.R.,  "The  Effectiveness  of  Base-Bleed  in  Reducing  Drag  of 
Boattailed  Bodies  at  Supersonic  Velocities,"  Memorandum  Report  No.  1244,  U.S. 
Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland 
21005,  1960  (AD  234315). 
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drag  reduction  has  been  studied.  Thus,  only  limited  attention  has  been  focus¬ 
sed  on  the  'base  bleed'  problem  in  transonic  flow.  A  limited  study  made  in 
the  transonic  flow  regime  has  been  reported  in  Reference  5  which  describes  the 
effects  of  base  bleed  on  various  afterbody  configurations. 

Most  of  the  work  using  the  'base  bleed'  concept  has  been  either  experi¬ 
mental  or  semi- empirical  in  nature.  Sophisticated  numerical  techniques  have 
not  yet  been  utilized  to  predict  the  effects  of  base  bleed  on  the  base  drag 
reduction.  Limited  computational  work  was  reported  recently  by  Sullins,  et 
al.6  Their  work  dealt  with  the  numerical  computation  of  the  base  region  flow 
of  a  supersonic  combustion  ramjet  engine  using  two-dimensional  Navier-Stokes 
equations.  They  computed  the  flow  field  in  the  vicinity  of  the  base  with 
parallel  gas  injection  and  established  the  effect  of  base  injection  on  such 
flows. 

Because  of  the  recent  advances  in  computer  technology,  numerical  computa¬ 
tional  capabilities  have  been  developed  to  predict  the  aerodynamic  behavior  of 
artillery  shell.  Recent  papers7*0  have  reported  the  development  and  applica¬ 
tion  of  the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  computational  techni¬ 
que  to  predict  the  flow  about  slender  bodies  of  revolution  at  transonic 
speeds.  This  technique  has  been  modified  for  base  flow  analysis  and  the 
resulting  new  numerical  capability9  is  used  here  to  predict  the  base  pressure 
of  shell  at  transonic  speeds  including  the  effect  of  base  bleed.  Computed 
results  show  quantitative  and  qualitative  details  of  the  base  flow  structure. 
The  technique  computes  the  full  flow  field  over  the  projectile  at  transonic 
speeds;  therefore,  all  three  components  of  the  total  drag  (pressure,  viscous, 
and  base  drag)  are  computed.  This  computational  technique  is  then  applied  to 
predict  the  effects  of  base  bleed  on  the  base  drag  reduction  at  transonic 


5.  Sykes ,  D.M.,  " Cylindrical  and  Boattailed  Afterbodies  in  Transonic  Flow 
with  Gas  Ejection,"  AIAA  Journal,  Vol .  8,  No.  3 ,  March  1970 ,  pp.  588-589 . 

6.  Sullins,  G.A.,  Anderson,  J.D.,  and  Drummond,  J.P.,  " Numerical  Investiga¬ 
tion  of  Supersonic  Base  Flow  with  Parallel  Injection,"  AIAA  Paper  No.  82- 
1001,  June  1982. 

7.  Nietubicz,  C.J.,  Pulliam,  T.H.,  and  Steger,  J.L.,  "Numerical  Solution  of 
the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  Equations,  "  ARBRL-TR- 
02227,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland  21005,  March  1980  (AD  A085716). 

'  8 .  Nietubicz,  C.J.,  "Navier-Stokes  Computations  for  Conventional  and  Hollow 
Projectile  Shapes  at  Transonic  Velocities,"  ARBRL-MR-03184,  U.S .  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland  21005, 
July  1982  (AD  A116866). 

9.  Sahu,  J.,  Nietubicz,  C.J.,  and  Steger,  J.L.,  "Numerical  Computation  of 

Base  Flow  for  a  Projectile  at  Transonic  Speeds,"  ARBRL-TR-02495,  U.S.  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland  21005, 
June  1983  (AD  A130293). 
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speeds.  The  combined  effect  of  boattailing  and  base  bleed  is,  however,  not 
considered  in  this  report. 

Brief  descriptions  of  the  physical  problem  and  the  governing  equations 
are  given  in  Sections  II  and  III,  respectively.  The  computational  technique 
and  the  method  of  solution  are  discussed  in  Section  IV.  In  Section  V,  results 
are  shown  for  transonic  base  pressure  computations  for  a  6-caliber  secant- 
ogive-cylinder  shape  for  .9  <  M  <  1.2,  with  and  without  base  bleed.  Velocity 
vector  plots  and  stream  function  contour  plots  are  presented  to  show  the  qual¬ 
itative  features  of  the  flow  field  in  the  base  region.  Quantitative  compari¬ 
sons  of  base  drag  and  the  total  drag,  both  with  and  without  base  injection, 
have  been  made.  The  encouraging  results  show  that  the  present  computational 
technique  can  be  used  to  study  the  effects  of  base  bleed  on  base  drag  and  thus 
can  have  a  positive  impact  on  the  XM864  devleopment.  Although  results  in  this 
report  are  presented  for  transonic  speeds,  current  computational  efforts  are 
directed  at  supersonic  velocities. 


II.  PHYSICAL  PROBLEM 

The  physical  problem  deals  with  the  transonic  flow  over  a  projectile, 
including  the  base  region.  Although  the  entire  projectile  flow  is  computed, 
the  emphasis  is  on  the  flow  field  in  the  base  region  of  the  projectile.  A 
small  amount  of  air  is  injected  at  the  projectile  base  in  the  direction  paral¬ 
lel  to  the  primary  flow.  The  injection  at  the  base  can  be  concentrated  at  the 
center  of  the  base  or  spread  throughout  the  entire  base.  In  the  present  work, 
however,  the  injection  takes  place  over  90%  of  the  base.  Figure  1  shows  a 
schematic  illustration  of  the  base  region  flow  field  with  base  injection.  The 
dividing  streamline  separates  the  recirculary  base  flow  from  the  primary 
external  flow.  The  flow  field  is  dominated  by  separation  and  mixed  regions  of 
locally  supersonic  and  subsonic  flows. 

The  complete  set  of  time-dependent  generalized  axisymmetric  thin-layer 
Navier-Stokes  equations  is  solved  to  obtain  a  numerical  solution  to  this 
problem.  The  numerical  technique  used  is  an  implicit  finite-difference 
scheme.  Although  time-dependent  calculations  are  made,  the  transient  flow  is 
not  of  primary  interest  at  the  present  time.  The  steady  flow  is  the  desired 
result  which  is  obtained  in  a  time  asymptotic  fashion. 


III.  GOVERNING  EQUATIONS 

The  Azimuthal  Invariant  (or  Generalized  Axisymmetric)  thin-layer  Navier- 
Stokes  equations  for  general  spatial  coordinates  5,  n,  ±  can  be  written  as7 


where 


3S 

3C 


K  =  5(x,y,z,t)  is  the  longitudinal  coordinate 
n  =  n(y,z,t)  is  the  circumferential  coordinate 


(1) 
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5  =  c(x,y,z,t)  is  the  near  normal  coordinate 
t  =  t  is  the  time 


The  vector  of  dependent  variables  q  and  the  flux  vectors  E,  G,  H  are 
given  as 


—  “ 

p 

pu 

pu 

a  i 

E  =  J  1 

puU  +  £xP 

PV 

pvU  +  £  p 

» 

y 

pw 

pwU  +  5zp 

e 

_ 

(e  +  p)  U  -  CtP 

pW 

puW  +  cxP 

PVW  + 

pwW  +  tzP 

(e  +  p)  W  -  ttp 
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*  -1 
H  =  J  x4> 

n 


0 

0 


pV[R^(u  -  et)  +  RC(W  -  ct)3 
-pVR^n(V  -  nt)  -  p/(R<t>n) 


where  J  is  the  Jacobian  of  transformation. 

The  thin  layer  viscous  terms  valid  for  high  Reynolds  number  flow  are 

A 

contained  in  the  vector  Sf  where 


S  = 


0 


"<‘x  +  S2  +  cz)uc  +  Wl'Vt  *  Vc  *  ‘z\>4x 

uUJ  +  C*  +  c|)vc  +  (u/3)(cxuc  +  cyv^  +Cz«c)cy 

pUx  +  Qy  +  f  (u/3)(cxuc  +  cyvc  +  czwc)cz 

(U2  +  C*  +  Cp[0.5u(u2  +  v2  +  w2)^  +  <Pr_1  (y  -  l)*1  (a2) 

+  (u/3)(cxu  +  c  v  +  czw)  x  (cxuc  +  cyvc  +  czwc)} 
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The  velocities 


u  -  «t  +  5xu  +  y  *  h* 

V  =  nt  +  nu  +  nv  +  n  w  (2) 

t  x  y  z 

W  -  {t  +  V  *  V  +  'zW 

represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  nondimensional ized  with 
respect  to  a^  (the  free  stream  speed  of  sound).  The  density  (p)  is  referenced 
to  and  total  energy  (e)  to  p^a^.  The  local  pressure  is  determined  using 
the  equation  of  state. 


p  *  (y  -  l)[e  -  0.5p(u2  +  v2  +  w2)] 


(3) 


where  y  is  the  ratio  of  specific  heats. 

In  Equation  (1)  a  thin-layer  approximation  is  used  and  the  restrictions 
for  axisymmetric  flow  are  imposed.  The  details  can  be  found  in  References  8 
and  9  and  are  not  discussed  here.  Equation  (1)  contains  only  two  spatial 
derivatives;  however,  it  retains  all  three  momentum  equations,  thus  allowing  a 
degree  of  generality  over  the  standard  axisymmetric  equations.  In  particular, 
the  circumferential  velocity  is  not  assumed  to  be  zero,  thus  allowing  computa¬ 
tions  for  spinning  projectiles  or  swirl  flow  to  be  accomplished.  There  is 
some  evidence  which  indicates  that  base  pressure  can  change  due  to  the  spin  of 
a  projectile.  Although  the  present  work  considers  base  flow  with  no  spin, 
base  flow  with  spin  is  of  interest  and  can  be  studied  using  the  present 
technique. 

For  the  computation  of  turbulent  flows  a  turbulence  model  must  be  sup¬ 
plied.  In  the  present  calculations  a  Cebeci-type  two  layer  algebraic  eddy 
viscosity  model  as  modified  by  Baldwin  and  Lomax10  is  used.  In  their  two 
layer  model  the  inner  region  follows  the  Prandtl-Van  Driest  formulation. 
Their  outer  formulation  can  be  used  in  wakes  as  well  as  in  attached  and  sepa¬ 
rated  boundary  layers.  In  both  the  inner  and  outer  formulations  the  distribu¬ 
tion  of  vorticity  is  used  to  determine  length  scales,  thereby  avoiding  the 
necessity  of  finding  the  outer  edge  of  the  boundary  layer  (or  wake).  The 
magnitude  of  the  local  vorticity  for  the  axisymmetric  formulation  is  given  by 


0) 


- '  O2  *  < 


3v 

dZ 


3W\ 2  +  ( dW  3U  \  2 
3y'  'ax  "  3z; 


(4) 


10.  Baldioin ,  B.S.,  and  Lomax,  H.t  "Thin-Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flou>8,"  AIAA  Paper  A Jo.  78-257 ,  1978 . 
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In  determining  the  outer  length  scale  a  function10 

F (y)  ■  y|“>|  Cl  -  exp(-y+/A+)] 


(5) 


is  used  where  y+  and  A+  are  the  conventional  boundary  layer  terms.  For  the 
base  flow  (or  wake  flow)  the  exponential  term  of  Equation  (5)  is  set  equal  to 
zero. 


IV.  NUMERICAL  METHOD 


A.  Computational  Algorithm 

An  implicit,  approximate  factorization,  finite-difference  scheme  in  delta 
form  as  described  by  Beam  and  Warming11  is  used.  An  implicit  method  was 
chosen  because  it  permits  a  time  step  much  greater  than  that  allowed  by  explic¬ 
it  schemes.  For  problems  in  which  the  transient  solution  is  of  no  interest, 
this  offers  the  advantage  of  being  able  to  reach  the  steady  state  solution 
faster  than  existing  explicit  schemes. 

The  Beam-Warming  implicit  algorithm  has  been  used  successfully  in  various 
applications.7”13  The  algorithm  can  be  first  or  second  order  accurate  in  time 
and  second  or  fourth  order  accurate  in  space.  The  equations  are  factored 
(spatially  split)  which  reduces  the  solution  process  to  one-dimensional  prob¬ 
lems  at  a  given  time  level.  Central  difference  operators  are  employed  and 
the  algorithm  produces  block  tridiagonal  systems  for  each  space  coordinate. 
The  main  computational  work  is  contained  in  the  solution  of  these  block  tri¬ 
diagonal  systems  of  equations. 


11.  Beam,  R.,  and  Warming,  R.F.,  "An  Implicit  Factored  Scheme  for  the 

Compressible  Navier-Stokes  Equations,"  AIAA  Journal,  Vol .  16,  No.  4, 
April  1978 ,  pp.  393-402.  ' 

12.  Stegev,  J.L.,  "Implicit  Finite  Difference  Simulation  of  Flow  About 

Arbitrary  Geometries  with  Application  to  Airfoils,"  AIAA  Journal,  Vol. 
16,  No.  7,  July  1978,  pp.  679-686.  ” 

13.  Pulliam,  T.H.,  and  Steger,  J.L.,  "On  Implicit  Finite-Difference  Simula¬ 
tions  of  Three-Dimensional  Flow,  "  AIAA  Journal,  Vol.  18,  No.  2,  February 
1980,  pp.  159-167. 
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B.  Finite  Difference  Equations 


The  resulting  finite  difference  equations,  written  in  delta  form  are 
(I  +  h6^An  -  Ejj'1?  A  J)(I  +  hS^C"  -  EjJ-1 

-  hRe’^J"1^)  (qn+1  -  qn)  =  -At(6^En  +  6^Gn  (6) 

-Re-1  fi^Sn)  -  AtHn  -  eEJ_1[(7^)2  +  (V^)2]Jqn 


Here  h  =  At  because  only  first  order  accuracy  in  the  time  differencing  is 
needed  for  the  steady  state  flows  considered  here.  This  choice  corresponds  to 
the  Euler  implicit  time  differencing.  The  6's  represent  central  difference 
operators,  a  and  7  are  forward  and  backward  difference  operators, 

A  A 

A 

respectively.  The  Jacobian  matrices  A  =  — ,  C  =  — *■  together  with  the  coeffi- 

.  .  .  3q  3q 

ciertt  matrix  M  obtained  from  the  local  time  linearization  of  S  are  described 
in  detail  in  Reference  6.  Fourth  order  explicit  (e^)  and  implicit  ( e j )  numer¬ 
ical  dissipation  terms  are  incorporated  into  the  differencing  scheme  to  damp 
high  frequency  growth  and  thus  to  control  the  nonlinear  instabilities.  A 
typical  range  for  the  smoothing  coefficients  is  =  (1  to  5)  At  with 
£j  =  3e£. 

C.  Flow  Field  Segmentation 

Figure  2  is  a  schematic  illustration  of  the  flow  field  segmentation  used 
to  compute  the  entire  projectile  flow  field  including  the  base  flow.  It  shows 
the  transformation  of  the  physical  domain  into  the  computational  domain  and 
the  details  of  the  flow  field  segmentation  procedure  in  both  domains. 

The  cross  hatched  region  represents  the  projectile.  The  line  BC  is  the 
projectile  base  and  the  region  ABCD  is  the  base  region  or  the  wake.  The  line 
AB  is  a  computational  cut  through  the  physical  wake  region  which  acts  as  a 
repetitive  boundary  in  the  computational  domain.  Implicit  integration  is 
carried  out  in  both  5  and  c  directions.  (See  Figure  2.)  Note  the  presence  of 
the  lines  BC  (the  base)  and  EF  (nose  axis)  in  the  computational  domain.  They 
both,  however,  act  as  boundaries  in  the  computational  domain  and  special  care 
must  be  taken  in  forming  the  block  tridiagonal  matrix  in  the  5  direction.  The 
details  are  presented  in  the  next  section. 

D.  Implementation  of  Boundary  Conditions 

1 .  Base  Flow  without  Base  Injection. 

The  no-slip  boundary  conditions  for  viscous  flow  are  enforced  by 

setting 
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U  =  V  *  W  =  0 


(7) 


on  the  projectile  surface  except  for  the  base.  At  the  projectile  base  the 
velocity  component  normal  to  the  base  is  set  to  zero,  i.e.,  U  =  0,  while  other 
flow  variables  are  set  equal  to  those  at  grid  point  next  to  the  base. 

Along  the  computational  cut  (line  AB),  the  flow  variables  above  and 
below  the  cut  were  simply  averaged  to  determine  the  boundary  conditions  on  the 
cut.  On  the  centerline  of  the  wake  region,  a  symmetry  condition  is  imposed: 


3U 

3Z 


0 


3 V 
3z 


0 


w  =  0 


If  ■ ° 


(8) 


Free  stream  conditions  are  used  at  the  outer  boundary.  First  order 
extrapolation  for  all  flow  variables  is  used  at  the  downstream  boundary  (lines 
AD  and  AG).  During  transient  calculations,  use  of  a  specified  outflow  pres¬ 
sure  can  give  rise  to  numerical  oscillations  in  the  base  region  flow  field. 
Eventually,  these  grow  and  swamp  the  solution.  This  difficulty  is  avoided  by 
simply  extrapolating  pressure  to  the  downstream  boundary  which  is  the  proce¬ 
dure  always  used  with  supersonic  outflow.  A  combination  of  extrapolation  and 
symmetry  is  used  on  the  nose  axis  (line  EF). 

As  a  result  of  the  flow  field  segmentation  procedure  described 
earlier,  the  block  tridiagonal  matrix  in  the  5  direction  has  elements  at  J  = 
JB,  JB+1  which  are  treated  as  internal  boundaries  in  the  computational  domain. 
(J  =  JB  represents  the  projectile  base  and  J  =  JB+1  is  the  nose  axis.)  The 
block  tridiagonal  matrix  in  the  5  direction  takes  the  following  form  (after 
setting  £j  =  0  to  simplify  the  illustration) 
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—  — 

n 

-  — 

n 

— 

I  a3 

Aq2 

rhs2 

-a2  \  h 

Aq3 

t 

RHS3 

"aJB-2  1  AJB 

1 

0  I  0 

AqJB 

= 

0 

0  10 

iqJB+l 

0 

~aJB+1  \  AJB+3 

~AJMAX-2  1 

-  M 

AqJMAX-l 

RHSJMAX-1 

L. 

At  A 

Here  A's  denote  the  quantity  A  and  I  is  a  5x5  identity  matrix.  Note  the 

appearance  of  two  uncoupled  block  tridiagonal s.  The  rows  at  JB  and  JB+1  are 
particularly  simple  because  boundary  conditions  are  updated  explicitly  at  the 
end  of  inversions.  These  changes  were  easily  implemented  in  a  modular  fashion 
into  an  existing  code  for  projectile  base  flow  computations.  One  simply  fills 
the  block  tridiagonal  matrix  ignoring  the  base  JB  and  the  nose  axis  JB+1. 
Elements  in  these  rows  are  then  overloaded  as  shown  above.  The  flow  field 
segmentation  does  not  affect  the  block  tridiagonal  matrix  in  the  t  direction. 

2.  Base  Flow  with  Base  Injection. 

The  boundary  conditions  used  for  base  flow  with  mass  addition  are 
presented  here.  The  boundary  conditions  along  the  projectile  surface,  at  the 
cut  and  downstream  boundary  all  remain  the  same  as  previously  described. 
Along  the  base  boundary  the  following  conditions  are  imposed: 


u  =  u. 


V  -  Vj  =  0 


w  =  wj  -  wjg_i  (grid  point  next  to  the  base) 

tkhu'.cm.  kpmts  ssmos 

STinfO  BRANCH^ 

BLOB.  305 


p  =  Pj  =  Pst 


The  stagnation  density  is  obtained  from  the  following  relation. 
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(10) 


-  (1  ♦  ^ 


The  amount  of  air  injected  Into  the  base  region  can  be  specified  by  the  mass 

flow  rate,  m..  Since  p.  and  A,-  are  known,  can  be  calculated  for  any  given 
J  J  J  J 

mass  flow  rate.  Rather  than  specifying  m., however,  it  is  customary  to  specify 

^  • 

a  mass  injection  parameter,  Ij  where  Ij  =  Most  of  the  results  with 

base  bleed  are  presented  in  terms  of  this  parameter  in  the  next  section.  It 
is  important  to  remember  that  the  smaller  the  mass  injection  parameter,  the 
smaller  is  the  amount  of  mass  injected  at  the  base. 

E.  Computational  Grid 


The  finite  difference  grid  used  for  the  numerical  computations  was 
obtained  from  a  grid  generator  developed  and  reported  in  Reference  14.  This 
program  allows  arbitrary  grid  point  clustering,  thus  enabling  grid  points  for 
the  projectile  shapes  to  be  clustered  in  the  vicinity  of  the  body  surface. 
The  grid  consists  of  108  points  in  the  longitudinal  direction  and  50  points  in 
the  radial  direction.  The  full  grid  is  shown  in  Figure  3  while  Figure  4  shows 
an  expanded  view  of  the  grid  in  the  vicinity  of  the  projectile.  The  computa¬ 
tional  domain  extended  to  4  body  lengths  in  front,  4  body  lengths  in  the 
radial  direction  and  4  body  lengths  behind  the  base  of  the  projectile.  The 
grid  points  in  the  normal  direction  were  exponentially  stretched  away  from  the 
surface  with  the  minimum  spacing  at  the  wall  of  .00002D.  This  spacing  locates 
at  least  two  points  within  the  laminar  sublayer. 

The  grid  shown  in  Figure  4  was  generated  in  two  segments.  First,  the 
grid  in  the  outer  region  is  obtained  using  an  elliptic  solver14  for  the  ogive 
portion  and  straight-line  rays  for  the  remaining  portion  which  runs  all  the 
way  to  downstream  boundary.  Second,  the  grid  in  the  base  region  is  obtained 
simply  by  extending  the  straight  lines  perpendicular  to  line  AB  down  to  the 
center  line  of  symmetry  (line  CO).  An  exponential  stretching  with  the  minimum 
spacing  of  .000020  at  line  AB  is  used.  It  should  be  noted  that  the  same  mini¬ 
mum  spacing  .000020  is  specified  on  both  sides  of  the  cut  thus  maintaining  a 
smooth  variation  of  grid  across  the  cut.  This  spacing  could,  of  course,  be 
increased  downstream  of  the  base.  The  number  of  grid  points  above  and  below 
line  AB  is  the  same  (50  points).  As  can  be  seen  in  Figure  4,  the  grid  points 
are  clustered  near  the  nose-cylinder  junction  and  at  the  projectile  base  where 
appreciable  changes  in  flow  variables  are  expected. 

As  indicated  in  Figure  4,  the  fine  viscous  grid  follows  the  cut  labeled 
as  AB  in  Figure  2.  Insofar  as  the  viscous  shear  layer  begins  to  neck-down 
shortly  behind  the  base,  much  of  this  fine  grid  resolution  is  wasted.  As  a 


14 .  Steger ,  J.L.,  Nietubicz ,  C.J,,  and  Heavey,  K.R.,  " A  General  Curvilinear 
Grid  Generation  Program  for  Projectile  Configurations,"  ARBRL-MR-03142 , 
U. 5.  Army  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
Maryland  21005 ,  October  1981  (AD  A1 07334). 
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consequence,  logic  has  been  implemented  to  adjust  the  grid  cut  AB  to  the 
viscous  shear  layer.  Such  a  grid  is  shown  In  Figure  5  in  which  the  height  of 
the  cut  is  determined  from  a  moment  of  shear  subject  to  various  constraints 

and  averaging.  Specifically,  the  cut  height,  7,  at  each  J-location  is  deter¬ 
mined  by  the  relation 


-  £(6z  UJL>2  ZJL  +  eD/2 

ZJ  =  - E(6z  'ujl7z  +"e - 


(11) 


where  the  summation  is  carried  out  only  for  those  points  within  an  interval 
.2D  <  Zj^  <  D/2.  Here  D  is  the  base  diameter,  &z  is  a  central  difference 

operator  and  e  is  a  positive  parameter  which  ensures  a  standard  grid  if  all 
<5z  Ujl  are  zero  or  if  e  is  very  large.  «.  Additional  averaging  is  used  in  the  x- 

direction  (longitudinal  direction).  Preliminary  results  have  been  obtained 
using  the  grid  shown  in  Figure  5  and  further  computations  are  underway. 


V.  RESULTS 

The  model  geometry  used  in  the  present  study  is  shown  in  Figure  6.  It 
consists  of  a  3  caliber  secant-ogive  nose  and  a  3  caliber  cylinder. 

The  free  stream  Reynolds  number  for  the  series  of  computations  was  fixed 
at  4.5  x  106  based  on  the  total  model  length.  The  computations  were  started 
from  free  stream  conditions  and  marched  in  time  to  obtain  the  steady  state 
solution.  The  initial  calculation  was  made  for  M  =  0.9.  Previously  converged 
solutions  were  then  used  as  starting  conditions  for  additional  Mach  number 
runs  to  achieve  faster  convergence.  The  results  are  now  presented  for  both 
(i)  base  flow  without  injection  and  (il)  base  flow  with  injection. 

Figures  7  and  8  show  the  distribution  of  the  surface  pressure  coeffi¬ 
cient,  Cp  as  a  function  of  axial  position  without  and  with  mass  injection  at 

the  base,  respectively.  The  value  of  Cp  beyond  X/D  =  6  is  the  value  of  pres¬ 
sure  coefficient  along  the  line  extending  from  the  cylinder  portion  straight 
to  the  downstream  boundary.  When  there  is  no  mass  injection  at  the  base  the 
pressure  distribution  in  Figure  7  reflects  the  shock  pattern  that  typically 
occurs  on  shell  at  transonic  velocities,  the  rapid  expansion  which  occurs  at 
the  blunt  base  and  the  recompression  that  occurs  downstream  of  the  base.  The 
pressure  coefficient  distribution  for  a  ease  with  large  mass  addition  is  shown 
in  Figure  8.  The  previously  observed  rapid  expansion  at  the  base  and  recom¬ 
pression  downstream  of  it  are  seen  to  be  virtually  eliminated. 

Figure  9  shows  the  velocity  vector  field  in  the  base  region  for  M  =  0.9, 
a  «  0  and  Ij  =  0.  Each  vector  shows  the  magnitude  and  the  direction  of  the 

velocity  at  that  point.  The  figure  shows  the  velocity  field  when  there  is  no 
base  bleed  and  the  reci rculatory  flow  in  the  base  region  is  clearly  evident. 
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The  velocity  vector  plots  in  Figures  10,  11  and  12  show  the  effect  of 
base  bleed  on  the  near  wake  flow  field.  Figure  10  shows  the  effect  of  base 
bleed  for  the  case  when  the  mass  injection  parameter  is  very  small  (Ij  =  .01). 

The  change  in  the  flow  field  is  not  very  dramatic.  In  Figure  11  the  mass 
injection  parameter  is  increased  to  .07  and  the  effect  of  base  bleed  can  be 
clearly  seen.  The  near  wake  flow  field  has  changed  considerably.  Figure  12 
shows  the  effect  of  base  bleed  for  a  higher  mass  injection  parameter,  Ij  = 

.13.  The  flow  field  in  the  base  region  has  now  been  dramatically  altered. 
The  recirculation  pattern  has  been  totally  eliminated. 

The  next  four  Figures  (13,  14,  15  and  16)  are  stream  function  contour 
plots  in  the  wake  region,  again  for  M  =  0.9  and  a  =  0.  All  these  figures  are 
deliberately  stretched  in  y  direction  (not  drawn  to  the  same  scale  in  x  and  y) 
to  show  the  flow  pattern  in  the  base  region  as  clearly  as  possible.  Figure  13 
is  for  the  case  of  base  flow  with  no  mass  injection.  It  clearly  shows  the 
recirculation  region  and  the  position  of  the  dividing  streamline  which  separ¬ 
ates  the  reci rcul atory  base  flow  from  the  main  flow.  The  reattachment  point 
is  about  2  calibers  down  from  the  base.  Note  the  strong  shear  layer  in  the 
base  region. 

Figures  14,  15  and  16  show  the  flow  pattern  in  the  base  region  with  mass 
injection  allowed  at  the  base.  Figure  14  shows  the  effect  of  base  bleed  when 
the  mass  injection  parameter  is  very  small  ( I j  =  .01).  The  reattachment  point 

remains  at  about  the  same  place  as  with  no  injection  at  the  base.  The  flow 
pattern  has  changed  slightly  as  can  be  seen  by  the  dividing  streamline;  how¬ 
ever,  the  recirculation  region  has  not  changed  dramatically.  In  Figure  15,  the 
mass  injection  parameter,  I j ,  has  been  increased  to  .07  and  now  the  effect  of 

mass  injection  can  be  clearly  seen.  The  reattachment  point  has  moved  further 
downstream.  The  flow  pattern  in  the  near  wake  flow  field  has  changed  consid¬ 
erably  and  the  separation  bubble  is  reduced  in  size.  When  the  mass  injection 
parameter  is  increased  further,  Ij  =  .13,  its  effect  on  the  flow  field  in  the 

base  region  is  apparent.  Figure  16  shows  that  dramatic  change  in  the  flow 
field.  The  recirculation  region  has  been  eliminated  and  the  shear  layer  has 
been  displaced  markedly.  Eliminating  the  recirculation  in  the  wake,  increases 
the  base  pressure  and  thus  reduces  the  base  drag. 

A  more  critical  look  at  the  computational  results  is  presented  in  Figures 
17  through  20.  These  figures  show  the  quantitative  details  of  projectile  flow 
field.  Figure  17  shows  the  variation  of  base  drag  with  mass  injection  rates 
for  M  =  0.9  and  a  =  0.  The  reduction  in  base  drag  with  base  injection  can  be 
clearly  seen.  The  reduction  in  base  drag  increases  with  the  increase  in  the 
injection  rate. 

Since  the  entire  projectile  flow  field,  including  the  base  flow,  has  been 
computed,  all  three  drag  components  have  been  computed  and  thus  the  total  drag 
is  determined.  Figure  18  shows  the  variation  of  the  total  drag  with  varying 
mass  injection  rates.  Again,  the  reduction  in  the  total  drag  is  apparent.  As 
the  injection  rate  is  increased,  the  reduction  in  total  drag  increases. 

Figures  19  and  20  show,  respecti vely ,  the  variation  of  base  drag  and  the 
total  drag  with  Mach  number  both  with  and  without  base  injection.  In  both  of 
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these  figures  the  computational  results  without  base  injection  are  shown  by 
the  solid  line  whereas  the  dotted  line  represent  the  computational  results 
obtained  with  injection.  The  reduction  in  base  drag  and  thus  total  drag  due 
to  base  injection  can  be  clearly  seen.  Figure  19  indicates  that  the  reduction 
in  base  drag  has  increased  with  an  increase  in  Mach  number  from  .9  to  .98 
while  from  M  =  1.0  to  1.2,  the  drag  reduction  is  apparently  constant.  In  both 
of  the  figures,  the  expected  drag  rise  in  the  transonic  speed  regime  is  well 
predicted  for  .9  <  M  <  1.2  and  the  reduction  in  base  drag  and  the  total  drag, 
due  to  base  bleed  has  been  clearly  demonstrated. 


VI.  SUMMARY 

A  promising  computational  capability  has  been  developed  which  computes 
the  full  projectile  flow  field,  including  the  reci rcul atory  base  flow,  at 
transonic  speeds  both  with  and  without  base  injection. 

Numerical  computations  have  been  made  for  Mach  numbers  .9  <  M  <  1.2  to 
predict  the  base  drag  and  the  total  drag  with  and  without  base  bleed.  Compu¬ 
ted  results  show  the  qualitative  features  of  the  flow  field  in  the  near  wake 
for  both  cases.  The  effect  of  base  injection  on  the  qualitative  nature  of 
base  flow  has  been  clearly  shown.  Quantitative  comparisons  of  base  drag  and 
the  total  drag  both  with  and  without  base  injection  have  been  made  with  each 
other.  For  M  =  0.9  and  a  =  0  the  computational  results  show  the  reduction  in 
base  drag  and  the  total  drag  for  several  mass  injection  parameters.  Results 
are  also  presented  for  .9  <  M  <  1.2  for  a  given  mass  injection  rate  and  the 
reduction  in  base  dr^g  and  the  total  drag  has  been  demonstrated  for  this  range 
of  transonic  speeds. 

Current  efforts  are  directed  at  the  numerical  computation  of  base  flow  at 
supersonic  speeds.  The  encouraging  results  obtained  thus  far  at  transonic 
speeds  indicate  that  the  computational  technique  shows  the  promise  of  predict¬ 
ing  the  base  drag  and  hence  the  total  drag  both  with  and  without  base  injec¬ 
tion.  Future  computational  efforts  will  investigate  the  combined  effect  of 
boattailing  and  base  bleed  on  the  total  aerodynamic  drag. 
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Figure  1.  Schematic  Illustration  of  Base  Region  Flow  Field  with  Base  Bleed 
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OUTFLOW 


gure  2.  Schematic  Illustration  of  Flow  Field  Segmentation 


OUTFLOW 


Figure  6.  Model  Geometry 
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Figure  7.  Longitudinal  Surface  Pressure  Distribution, 
M  =  0.9,  a  =  0,  Ij  =  0  (without  Base  Bleed) 
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Figure  8.  Longitudinal  Surface  Pressure  Distribution, 
M  =  0.9,  a  =  0,  Ij  =  .13  (with  Base  Bleed) 
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Figure  9.  Velocity  Vector  Field,  M  =  0.9,  a  =  0,  I j  =  0 
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Figure  13.  Stream 
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Figure  19.  Variation  of  Base  Drag  Coefficient  with  Mach  Number, 

a  =  0  (with  and  without  Base  Bleed) 
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Figure  20.  Variation  of  Total  Drag  Coefficient  with  Mach  Number, 

a  =  0  (with  and  without  Base  Bleed) 
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LIST  OF  SYMBOLS 


a 

a„ 

00 

A 

Aj 

CDt 


A  A 


Dt 
e 

E,  F,  q 

A 

H 

I 

*3 

J 

m. 

J 

M 

Mm 

p 


Pr 

R 

Re 

A 

s 

t 


u,v,w 


speed  of  sound 
free  stream  speed  of  sound 
cross-sectional  area  at  the  base 
injection  area  for  base  bleed 
base  drag  coefficient,  2  D^p^u^A 
specific  heat  at  constant  pressure 
pressure  coefficient,  2(p  -  pm)/M^ 
body  diameter  (57.15rmn) 
base  drag 

total  energy  per  unit  volume/p^a^ 

flux  vector  of  transformed  Navier-Stokes  equations 

n-invariant  source  vector 

identity  matrix 

mass  injection  parameter, 

\J 

Jacobian  of  transformation 

mass  flow  rate  for  air  injection  at  the  base,  p-u-A- 

J  «J 

Mach  number 

free  stream  Mach  number 
pressure/p^ 
free  stream  pressure 
Prandtl  number,  C/k 

oo  rj  oo 

body  radius 

Reynolds  number,  P00^0P/vm 
viscous  flux  vector 
physical  time 

Cartesian  velocity  component s/a^ 
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utt  free  stream  velocity 

U,V,W  Contravarlant  velocity  components/a^ 

x,y,z  physical  Cartesian  coordinates 

a  angle  of  attack 

y  ratio  of  specific  heats 

<  coefficient  of  thermal  conductivity/^ 

<  coefficient  of  thermal  conductivity  at  free  stream  conditions 

p  coefficient  of  viscosity/p^ 

u  coefficient  of  viscosity  at  free  stream  conditions 

oo 

€,n,c  transformed  coordinates  In  axial,  circumferential  and  radial 
di recti  on s 

p  density/p^ 

p  free  stream  density 

t  transformed  time 

$  circumferential  angle 

<S  central  difference  operator 

a  forward  difference  operator 

v  backward  difference  operator 

Superscript 

*  critical  value 


Subscript 
b  base 

j  jet  conditions 

J  longitudinal  direction 

L  normal  direction 

o  total  conditions 
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st 


stagnation  conditions 
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